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Abstract 

Joint distribution function of N eigenvalues of U (N) invariant random- 
matrix ensemble can be interpreted as a probability density to find N 
fictitious non-interacting fermions to be confined in a one-dimensional 
space. Within this picture a general formalism is developed to study the 
eigenvalue correlations in non-Gaussian ensembles of large random matri- 
ces possessing non-monotonic, log-singular level confinement. An effective 
one-particle Schrodinger equation for wave-functions of fictitious fermions 
is derived. It is shown that eigenvalue correlations are completely de- 
termined by the Dyson's density of states and by the parameter of the 
logarithmic singularity. Closed analytical expressions for the two-point 
kernel in the origin, bulk, and soft-edge scaling limits are deduced in a 
unified way, and novel universal correlations are predicted near the end 
point of the single spectrum support. 
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1 Introduction and basic relations 



Random matrices are the field-theoretical models which study the universal fea- 
tures of physical phenomena resulting from the symmetry constraints only. This 
is the reason why quite different physical problems get a unified mathematical 
description within the framework of the random-matrix theory |Q. In partic- 
ular, the applicability of so-called invariant matrix model given by the joint 
distribution function 

N N 

p N ({X})d{X} = — l[dX k \X k \ aP exp{-0v(X k )} J] |Ai-A/ (1) 

N k=l i>j=l 

of N eigenvalues {A} of a large NxN random matrix H ranges from the problem 
of two-dimensional gravity 0, through the spectral properties of the Dirac 
operator in quantum chromodynamics Q , to the mesoscopic electron transport 
in normal and superconducting quantum dots B. Here the eigenvalues {A} 
belong to entire real axis, — oo < {A} < +oo, and the partition function Zjq is 
determined from the normalization condition J pn = 1. Parameter 

(3 in Eq. ([l]) accounts for the symmetry of the problem, a is a free parameter 
associated with a logarithmic singularity, while v (A) is a non-singular part of 
the confinement potential 

V(A) = «(A)-alog|A|. (2) 

It is implied that dimension N of the matrix H is large enough, N 3> I. In this 
thermodynamic limit the matrix model Eq. (Q) becomes exactly solvable. 

Different physics, that is behind the model introduced, deals with different 
regions of spectrum that can be explored in the corresponding scaling limits. 
Up to now, the most study received the random-matrix ensemble Eq. ([!]) with 
U (N) symmetry ((3 = 2), where three types of universal correlations have been 
established in the origin g §, bulk §, H jl(J ||, and soft-edge gf, || scaling 
limits. Corresponding eigenvalue correlations are described by the universal 
Bessel, sine, and G— multicritical kernels, respectively. Various scaling limits of 
the model Eq. ([j]) have been investigated by using different methods, so that 
a unified treatment of the problem of spectral correlations in U (N) invariant 
ensembles is still absent. The purpose of this paper is to present a unified 
approach allowing us to explore the spectral properties of the U (N) invariant 
matrix model Eq. (|l]) with effective log-singular level confinement in an arbitrary 
spectrum range. 

The following representation [0 of the joint distribution function pn ({A}) 
is well-known in the random-matrix theory: 

PA r({A}) = |vI/ (A 1 ,...,A N )| 2 , (3) 
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*o (Ai, X N ) = — L det (Ai)||^ =1 ...jv ■ ( 4 ) 



'N\ 

As far as takes the form of the Slater determinant, pjv ({A}) can be thought 
of as a probability density to find N non-interacting fictitious fcrmions in the 
quantum states ipo, <fw-i a t the "spatial" points Ai,...,Ajv- The "wave- 
functions" of such fermions are uniquely determined by the set of polyno- 
mials P n (A) orthogonal on the entire real axis with respect to the measure 
dfx(X) = exp{-2F(A)}dA, 

Vn(A) =P„(A)exp{-^(A)} (5) 

so that the orthogonality relation 

dXip n (A) ip m (A) = 8 nm (6) 

holds. It follows from Eq. (Q) that the joint distribution function ({A}) can 
be represented as 

PN ({A}) - i det \\K N (X u A,)|| . , (7) 

where Km (A, A') (referred to as the "two-point kernel") 

K N (X,X') = J2<Pk(\)Vk(X')- (8) 

fe=0 

is completely determined by the wave- functions ip n . Due to an additional 
constraint on the wave-functions of three successive quantum states that re- 
sults from the recurrence equation Eq. ([l0|) below, only the highly excited 
states, ifN-i and ipw, contribute to the two-point kernel in accordance with the 
Christoffel-Darboux theorem |l4| : 

v ,m (AQ (fiN~i (A) - ip N (A) <pn-i (AQ 

AAr(A,A) = c7v A' — A ' ^ 

This formula simplifies significantly the mathematical calculations in the ther- 
modynamic limit N ^> 1. Effective Schrodinger equation for ipx, that is the 
cornerstone of our unified approach, will be derived in the next Section. 



2 Effective Schrodinger equation 

In the particular case of the Gaussian unitary ensemble (GUE) the wave- functions 
tp n (A) are well-known. They are eigenfunctions of a fermion confined by a 
parabolic potential fm. For general non-Gaussian ensemble Eq. ([I]) the cal- 
culation of such effective wave-functions can be done by an extension of the 
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Shohat's method ]l5|, that previously has been used by the authors Q to 
treat the problem of eigenvalue correlations in random-matrix ensembles with 
non-singular, strong level confinement. This method allows us to map a three- 
term recurrence equation 

AP„_i (A) = c n P n (A) + c„_ 1 P n _ 2 (A) (10) 

for polynomials P n (A) orthogonal on the entire real axis with respect to the 
measure dfx (A) = exp {—2V (A)} dX, 

+00 

(A) P n (A) P m (A) = S nm , (11) 

-oc 

onto a second-order differential equation for corresponding fictitious wave-functions 
ip„ . Coefficients c n appearing in Eq. ( |l0| ) are uniquely determined by the mea- 
sure dfi. 

In order to derive an effective Schrodinger equation, we note the following 
identity 

= A n (A) P„_! (A) - B n (A) P n (A) , (12) 

with functions A n (A) and B n (A) to be determined from the following consid- 
eration. Since dP n (A) /d\ is a polynomial of the degree n — 1, it can be repre- 
sented Jl4| through the Fourier expansion in the terms of the kernel Q n (t, A) = 
EtJ Pk (A) Pk (1) *S: 

*p»m _ Mt) *M) wy (13) 



dX J_ OD dt 
Integrating by parts in the last equation we get that 
dP n (A) 



dX 



f + °° ( dV dV\ 

2 J dn (t) Q„ (t, A) ^— - —J P n (t) . (14) 



Now, making use of the Christoffel-Darboux theorem, we conclude that unknown 
functions A n (A) and B n (A) in Eq. ( |l2|) are 

B n (A) = 2c„ £j ^2 - f ) P„ (,) ( t ) . (16) 

We also notice the identity that directly follows from Eqs. (|l^), (jl^), (jn]) and 
from oddness of dV/dX: 

B n (A) + B n -i (A) - -A n ^ (A) = -2^. (17) 
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Differentiating Eq. (|l2|), making use of the recurrence equation Eq. (|l0|), 
and bearing in mind relation Eq. (^) between P n (A) and tp n (A), one can obtain 
an exact differential equation for the wave- functions of fictitious fcrmions, that 
is valid for arbitrary n: 

*<P»M F n{ x)^ + g n (X)<p n (X) = o, (is) 



dX 2 " w dX 



where 



S„(A) = ^i + J^A„A„_ 1 -bJb„ + 2^ + ±^) (20) 
dX c„_i \ dX A n dX ) 

d 2 V fdV\ 2 1 dA n dV 



dX 2 \dXJ A n dX dX ■ 

Previously, equation of this type was known in the context of the random- 
matrix theory only for GUE, where V (A) = A 2 /2. For such a confinement 
potential both functions A n and B n can easily be computed from Eqs. ( |l5| ) and 
(fl6|), and are given by A„ (A) = 2c„ and B n (A) = 0. Taking into account that 
for GUE c„ = we end up with T n (A) = and Q n (A) = 2n + 1 - A 2 . This 

allows us to interpret tp n (A) as a wave-function of the fermion confined by a 
parabolic potential: 

\ 2 { ) + (2n + 1 - A 2 ) ^ UE (A) = 0. (21) 



In principle, the effective Schrodinger equation Eq. (|18j) applies to general non- 
Gaussian random-matrix ensembles as well, although the explicit calculation of 
T n (A) and Q n (A) in this situation may be a rather complicated task. However, 
significant simplifications arise in the thermodynamic limit n = N 1. 

To proceed with derivation of the asymptotic Schrodinger equation, we have 
to specify the form of confinement potential V introduced by Eq. (|^) . Choosing 
the regular part v (A) to be an even function, we set 

W Q )(A) = ^|A 2 ^alog|A|, (22) 
fc=i 

with dp > 0. The signs of the rest dk's can be arbitrary, allowing for non- 
monotonic level confining, but they should lead to an eigenvalue density sup- 
ported on a single connected interval (— -Dj\r, +Dn). Confinement potential 
y( Q ) (A) determines its own set of orthogonal polynomials P„ (A), and func- 
tions and Bn which are needed to construct an asymptotic second-order 
differential equation for the function ip^ (A) = | A| Q -Py (A) exp {— v (A)}. [Here 



5 



upper index a reflects the presence of the log-singular component in yw (A), 
and the restriction a > — \ takes place due to normalization Eq. (pT 



In accordance with Eqs. ( |15| ) and (^2|) it is convenient to represent A^' in 
the form 

A%^\)=AW(\) + aAW(\), (23) 

where 

«(A) = 2c„£"*iW(^»«)) 2 . (25) 
Analogously, Eq. (|l^) leads to the similar representation 

< ) (A)=S r W(A) + aSS(A), (26) 

with 

<iw=^r°°^<w^.w- (28> 



— oo 



In the above formulas and Sfig result from the regular component of 

confinement potential, while A^g an( i ^sing are caused by its log-singular part. 

First, it is easy to see that A^2 g (A) = due to evenness of the measure d\i. 
Second, the exact expression for immediately follows from the recurrence 



equation Eq. (10), whence we get 

Bffl(A) = i^£. (29) 

Calculation of the regular parts Ai^g and -Breg can be done along the lines 
presented in Ref. |l]|. Then, we immediately obtain that Affl (A) is expressed 
in terms of Dyson's density vd (A) as follows: 

^ (A) = * VD[X) , (30) 
^1 - (X/D N ) 2 



v (\\- 2 V & dv f31l 

where Dn = 2c m should be identified with the soft edge of the spectrum. [It 
is easy to see that a log-singular part of the confinement potential does not 



G 



contribute to the Dyson's density, so that in the thermodynamic limit there 
are no changes in Dn due to logarithmic singularity of confinement potential] . 
Expression for 5reg can be obtained by the use of the large— N version of the 
identity Eq. (0), that yields 



<» (A) = ^-A 

UN 



(a) m dv 



(32) 



whence 



(A) 



A TTV D (A) 



D 



N 



yji - (\/D N y 



dv 
dX 

(a) 



l-(-l) 



Now, having asymptotic representations for A y ^' and given by Eqs. fl30| ) 
and (|||), and taking into account Eqs. (fL8|), fll9| ) and (p0|), it is straightforward 
to obtain the following remarkable effective asymptotic Schrodingcr equation for 
the wave-functions tpff (A) = |A|" (A) exp {—v (A)} of highly excited states 
(N ^> 1) of fictitious fermions: 



(33) 



d 2 ^ ] 
dX 2 



+ 



d i 



~KV D (A) 



^l-(X/D N f 
7T I/jj (A) H — 



A ( Q ) 

dA 



^ (A) = 



(34) 



Also, due to Eq. ([12]), one can verify that the wave- functions of two successive 
quantum states are connected by the relationship 



7 (a) 

dA 



7TZ/23 (A) 



^1 - (A/Z^) 2 



(a) 



r(A)- 



_A_ 



/V 



^ (A) 



+ (- 



-1) (A)- 



(35) 



Equations ( p4|) and (|35j) provide a general basis for the study of eigenvalue 
correlations in non-Gaussian random-matrix ensembles in an arbitrary spectral 
range. In particular case of GUE, the Dyson's density of states is the celebrated 
semicircle, ^ UE (A) = n^ 1 \J D 2 N — A 2 with Dn = \/2N. The square-root law 



for i^5 UK (A) immediately removes the first derivative dipffi /dX in Eq. (|34|), 
providing us the possibility to interpret the fictitious fermions as those confined 
by a quadratic potential (q = 0). As far as the semicircle is a distinctive feature 
of density of states in GUE only, one will always obtain a first derivative in 
the effective Schrodinger equation for the non-Gaussian unitary ensembles of 
random matrices. Therefore, fictitious non-interacting fermions associated with 
non-Gaussian ensembles of random matrices live in a non-Hermitian quantum 
mechanics. 
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An interesting property of these equations is that they do not contain the 
regular part of confinement potential explicitly, but only involve the Dyson's 
density vd (analytically continued on the entire real axis) and the spectrum 
end point D^. In contrast, the logarithmic singularity (that does not affect the 
Dyson's density) introduces additional singular terms into Eqs. ( p4| ) and (|35|), 
changing significantly the behavior of the wave-function tpff near the origin 
A = 0. The influence of the singularity decreases rather rapidly outward from 
the origin. 

Structure of the effective Schrodinger equation leads us to the following fun- 
damental statements: (i) Eigenvalue correlations are stable with respect to non- 
singular deformations of the confinement potential, (ii) In the random-matrix 
ensembles with well-behaved confinement potential the knowledge of Dyson's 
density (that is rather crude one-point characteristics coinciding with the real 
density of states only in the spectrum bulk) is sufficient to determine the gen- 
uine density of states, as well as the n— point correlation function, everywhere. 
The latter conclusion is rather unexpected since it considerably reduces the 
knowledge required for computing n— point correlators. 



3 Local eigenvalue correlations 



Effective Schrodinger equation obtained in the preceding Section allows us to 
examine in a unified way the eigenvalue correlations in non-Gaussian ensembles 
with U (N) symmetry in different scaling limits. As we show below, it inevitably 
leads to the universal Bessel correlations in the origin scaling limit j^, |?J , to the 
universal sine correlations in the bulk scaling limit |8[ p| [l0| , |i"Tf , and to the 
universal G— correlations in the soft-edge scaling limit 13 . Corresponding two- 
point kernels are given by Eqs. (fio]), (ji^) and ([52]), respectively. 



3.1 Origin scaling limit and the universal Bessel law 

Origin scaling limit deals with the region of spectrum close to A = where 
confinement potential displays the logarithmic singularity. In the vicinity of 
the origin the Dyson's density can be taken as being approximately a constant, 
v d (0) = 1/Ajv (0), where A^r (0) is the mean level spacing at the origin in the 
absence of the logarithmic deformation of potential v. In the framework of this 
approximation, Eq. ( |34| ) reads 



d\ 2 \A%(0) A 2 



<p%> (A) - 0. (36) 



Solution to this equation that remains finite at A = can be expressed by means 
of Bessel functions: 
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fw+l ( A ) 



bV\J. 



7rA 

A(0) 
ttA 
A(0) 



(37) 
(38) 



where a and b are constants to be determined later, and A (0) — A2N (0) ~ 
A27V+1 (0). In accordance with Eq. (^), the two-point kernel can be written 
down as 



X AA' 



A' 



A 



A(0) 



ttA 
A(0) 



TTA 

A(0) 



ttA' 
A(0) 



(39) 



where the unknown factor c can be found from the requirement irfjv (A, A) = 
1/A(0). This immediately yields us c = — 7r/A(0). Defining now the scaled 
variable s — X s / A (0), we obtain that in the origin scaling limit the two-point 



kernel K olig (s, s') = limAr^oo 
law 



K { £ ] {\ s ,\ s ,)d\ s /ds 



takes the universal Bessel 



^orig (S, S J = - VSS' . (40) 



Formula (^0|) is valid for arbitrary a > — 4. Note, that a recent proof of 
universality of the Bessel kernel given in Ref. 7 was based on the Christoffcl 
theorem that imposed an artificial restriction on parameter a to be only 
positive integer. 



3.2 Bulk scaling limit and the universal sine law 

Bulk scaling limit has been explored in a number of works || |[ |l(| p| . It 
is associated with a spectrum range where the confinement potential is well 
behaved (that is far from the logarithmic singularity A = 0), and where the 
density of states can be taken as being approximately a constant on the scale 
of a few levels. In accordance with this definition one has 



-Kbuik (s, s') = lim K oiig (s, s') . 



s.s' — >oo 



(41) 



where s and s' should remain far enough from the end point Dn of the spectrum 
support. 

Taking this limit in Eq. (|40|), we arrive at the universal sine law 



i^buik (s, s') 



sin [ir (s — s') 
7r (s — s') 



(42) 



3.3 Soft-edge scaling limit and the universal G— multicritical 
law 

Soft-edge scaling limit is relevant to the tail of eigenvalue support where crossover 
occurs from a non-zero density of states to a vanishing one. It is known |l7| that 
by tuning coefficients dk which enter the regular part v of confinement potential 
[see Eq. (gg)], one can obtain a bulk (Dyson's) density of states which possesses 
a singularity of the type 

^2 \ m+l/2 / ^ 



^^l 1 -^ (43) 

with the multicritical index m = 0,2,4, etc., and TZn being a well-behaved 
function with 1Z n (±1) 7^ 0. [Odd indices m are inconsistent with our choice that 
the leading coefficient d p , entering the regular component v (A) of confinement 
potential, be positive in order to keep a convergence of integral for partition 
function Z N in Eq. (|)]. Such an 77i— th multicriticality can be achieved by 
many means, and the corresponding plethora of multicritical potentials V 1 -" 1 ^ is 
given by the equation 

dV (m *> (A) ^ [ +Dn dt / t 2 \ m+1/ \ ( t , 



So-called minimal multicritical potentials which correspond to 1Z n = const can 
be found in Refs. @ 0. 

Below we demonstrate that as long as multicriticality of order m is reached, 
the eigenvalue correlations in the vicinity of the soft edge become universal, 
and are independent of the particular potential chosen. The order m of the 
multicriticality is the only parameter which governs spectral correlations in the 
soft-edge scaling limit. 

Let us move the spectrum origin to its endpoint Dn, making the replacement 



D 



N 



1 + S ' 2 \nD N K N (l)) 



(45) 



that defines the m— th soft-edge scaling limit provided s <C (DnTZn (l)) 1 ^" oc 
TV 1 / 1 '*, with 

v*=m+-. (46) 

It is straightforward to show from Eqs. ( p4] ) and ( p5| ) that in the vicinity of 
the end point Djy the function (p^ (s) = tp£ (X — Djsr) obeys the universal 
differential equation 

Pit (*) - ^~^ ?n («) " * a(w "- 1) VJV (*) = 0, (47) 
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and that the following relation takes place: 

r-iV-* ( 2 



(fN-l (s) = <fN (s) + (-1) 



\TrD N n N (i) 



(p' N (s) . 



(48) 



Solution to Eq. ( [i"7|) which decreases at s — » +oo (that is at far tails of the 
density of states) can be represented through the function 



G(s|i/*) 



1 



-1/2 



7 -l(i~^) (^)^(i-*) 



[where -/^i^ l\ and 7_|_i a j \ are the Bessel functions] as follows: 

(p N ( s ) = aG (s\u*) , 



(49) 

s > 0, 
s < 0, 



(50) 



where a is an unknown constant. Making use of Eq. (48), we obtain that in the 
vicinity of the soft edge the two-point kernel is 

G (s|^*) G (sV) ■ si-' - G (gV) g (s|^*) ■ (sQ^ 

Aat (A s , A s /) = 6 ; , (51) 

s — s' 

where b is an unknown constant again. It can be found by fitting Jl2[ the density 
of states Kjsr (A s , A s ), Eq. (|5l|), to the Dyson's density of states z<d (A s ), Eq. 
(p3|), near the soft edge provided 1 <C s <C N 1 ^* . This yields us the value & = 
c^y 1 (ttcnTZn (l)) 1 ^" • Thus, we obtain that in the m— th soft-edge scaling limit, 
Eq. (f45|), the two-point kernel (s, s') = liniAr^oo [Kn (A s , A s /) dA^/ds] sat- 

isfies the universal law 



G(s|i/*)G' (s'|i/*)-sJ 



G(s'|i/*)G' (sK)-(s') 5 ~"* 



(52) 



These G— multicritical correlations are universal in the sense that they do 
not depend on the details of confinement potential, but only involve such an "in- 
tegral" characteristic of level confinement as the index m of the multicriticality. 
In particular case of m — 0, that is inherent in random-matrix ensembles with 
monotonic confinement potential, the function G coincides with the Airy func- 
tion, G (s||) = Ai (s), and the previously supposed universal Airy correlations 

PI 

Ai (s) Ai' (s') - Ai (a') Ai' (s) 



(53) 



are recovered. 
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It follows from Eq. that the density of states in the same scaling limit 



«/ffi(*)=(^G( B |i,*)) a -[G(sK)]V*-i (54) 

is also universal. The large— \s\ behavior of can be deduced from the known 
asymptotic expansions of the Bessel functions: 



dii + t^rlzi cos (Ml) 

7T + 4Tr\s\ \ V* J > 



s — > — oo, 



"£g00=< (55) 

Note that the leading order behavior as s — » — oo is consistent with the |s|" _1 
singularity of the bulk density of states, Eq. (|^J). 



4 Concluding remarks 

We have presented a general formalism for a treatment of the problem of eigen- 
value correlations in spectra of U (N) invariant ensembles of large random ma- 
trices with log-singular level confinement. An important ingredient of our anal- 
ysis is an effective one-particle Schrodinger equation [see Eqs. ( |l8| ) and (pi) )] 
for fictitious non-interacting fcrmions naturally appearing in the determinantal 
representation of the joint distribution function of N eigenvalues of large N x N 
Hermitian random matrix. The structure of the asymptotic equation Eq. ( p4[ ) 
allowed us to conclude that: (i) Eigenvalue correlations are stable with respect to 
non-singular deformations of confinement potential, (ii) In the random-matrix 
ensembles with well-behaved confinement potential the knowledge of Dyson's 
density (that is rather crude one-point characteristics coinciding with the real 
density of states only in the spectrum bulk) is sufficient to determine the gen- 
uine density of states, as well as the n— point correlation function, everywhere. 
We have also demonstrated that effective Schrodinger equation contains all the 
information about eigenvalue correlations in arbitrary spectrum range: the uni- 
versal Bessel kernel Eq. (fjo|) was found to describe eigenvalue correlations in the 
origin scaling limit; the universal sine kernel Eq. ( [4 4 ) was revealed in the bulk 
scaling limit; finally, we have shown that the soft-edge scaling limit is described 
by the novel universal G— multicritical kernel Eq. (|52|). 
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